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Abstract 

We consider the inverse problem of the detection of a single body, 
immersed in a bounded container filled with a fluid which obeys the sta- 
tionary Navier-Stokes equations, from a single measurement of force and 
velocity on a portion of the boundary. We obtain an estimate of stability 
of log- log type. 
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1 Introduction. 

In this paper we deal with an inverse problem associated to the stationary 
Navier-Stokes system. We consider a bounded set f2 C R" (we assume n = 
2,3, which are indeed the physically relevant cases) with a sufficiently smooth 
boundary <9fl This body is filled with a Navier-Stokes fluid of constant viscosity 
[i. We want to detect an object D immersed in this container, by collecting 
measurements of the velocity of the fluid motion and of the boundary forces, 
but we only have access to a portion T of the boundary dfl. Once a suitable 
boundary condition g is assigned on T, the velocity u — (u\, . . . ,u n ) and the 
pressure p of the fluid will obey the following Navier-Stokes system in tt\D: 

div a (u,p) = (u ■ V)u in Cl\D, 

div u =0 in 0\Z?, (11) 
u = g on r, 



■it 



= on dD. 



Here, 

a(u,p) = /i(Vu + Vu T ) — pi 
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is the stress tensor, where I denotes the n x n identity matrix. The ith com- 
ponent of the nonlinear term (u ■ V)u is given by 

3=1 3 

The last equation in (jl.ljl is called "no-slip condition" . Call v the outer normal 
vector held to dfl. Once g 6 H 3 (r) is assigned, one can, in principle, measure 
on r the induced normal component of the stress tensor 

%p = a(u,p) ■ v, (1.3) 

and try to recover D from a single pair of Cauchy data (g, ip) known on the 
accessible part of the boundary T. Under some additional regularity hypotheses, 
namely of <9£1 being of Lipschitz class, and g e Ha (r), the uniqueness for this 
inverse problem has been shown to hold (see [6]) by means of unique continuation 
techniques. This means that if ui and ui are two solutions of corresponding 
to an assigned boundary data g, for D = Di and D = D% respectively, and 
a(u%,pi) ■ v = o-(u2,P2) • v on T, then D\ — This paper is devoted to the 
analysis of the problem of stability, which we may roughly state as follows: 
Given two solutions (ui,pi) to (11.11) for two different Di, for i = 1,2, with the 
same boundary data g, if 

\\a(ui,pi) ■ v - a{u 2 ,V2) • HI < e > 
what is the rate of convergence of du(Di, D2) as e — » 0? 

(We denote by du the Hausdorff distance). In [7] we proved a rate of con- 
vergence of log-log type for the analogous problem in the simpler context of 
the Stokes system. Here, we will prove an equivalent result for the stationary 
Navier-Stokes equations. As for the (yet easier) Stokes problem, even if we add 
some a priori information on the regularity of the unknown domain, we can only 
obtain a weak rate of stability. This does not come unexpected since, even for 
much simpler problems of the same kind, the dependence of D from the Cauchy 
data is at most of logarithmic type. See, for example, [2] for a similar problem 
on electric conductivity, or |16j , |17j for an inverse problem regarding elasticity. 
There are, in fact, several counterexamples showing that the optimal rate of 
convergence for the inverse conductivity problem is no better than of log type 
(see p] and 0). 

The purpose of this paper is thus to prove a log-log type stability for the Haus- 
dorff distance between the boundaries of the inclusions, assuming a C 2,a regu- 
larity bound. Such estimates have been estabilished for various kinds of elliptic 
equations, for example, [2], [I], for the electric conductivity equation, [TB], [T7] 
for the elasticity system and the detection of cavities or rigid inclusions, and 
[7] for the Stokes equation. The main tool used to prove stability here and in 
the aforementioned papers [3J, [16], [17], [7] is a quantitative estimate of con- 
tinuation from boundary data, in the interior, in the form of a three spheres 
inequality (see Theorem 14. 1[) and its main consequences. However, while in [2j 
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the estimates are of log type for a scalar equation, here, and in [IB], J7] and 
0, only an estimate of log-log type could be obtained for a system of equations. 
To improve that, one would need a doubling inequality at the boundary for sys- 
tems of equations, which basically would allow to extend the reach of the unique 
continuation property up to the boundary. Unfortunately, to the present time, 
none are available; on the other hand they are known to hold in the scalar case. 
A very recent paper by Lin, Uhlmann and Wang (|14j) extended the validity 
of the three spheres inequality to linearized Navier-Stokes systems: this allows 
us to apply it to differences of solutions of see Proposition 14.51 In order 

to adapt this result to the Navier-Stokes equations, however, we are forced to 
add yet more a priori hypotheses on the solutions, mainly because of the non- 
linear character of the equations. Proposition 14.51 in fact, applies to linearized 
stationary Navier-Stokes systems with coefficients bounded in an appropriate 
norm. We meet this request by restricting the choice of boundary data, i.e. we 
require a strong regularity bound- i.e., C 1,Q -on the boundary data. 
The present paper has the same structure as [7] and closely follows [IB], [17] . 
The main ingredients are the following: 

1. An estimate of propagation of smallness from the interior. The proof of 
this estimate relies essentially on the three spheres inequality. Once we 
adapt it to , the estimate itself is adapted effortlessly. 

2. A stability estimate of continuation from the Cauchy data. This result 
also relies heavily on the three spheres inequality, but in order to obtain 
a useful estimate of continuation near the boundary, we need to extend a 
given solution of (jl.lj) to a larger domain, so that the extended solution 
solves a similar system of equations. We may then apply the stability 
estimates from the interior to the extended solution and treat them like 
estimates near the boundary for the original one. 

In Section 2, we state the apriori hypotheses we will need throughout the paper, 
recall some useful properties of the direct problem, and give the main result in 
Theorem 12.41 Section 3 contains the estimates of continuation in Propositions 
I3.ll and f3.21 and Propositions 13 . 31 and 13.41 These deal, in turn, with the stability 
estimates of continuation from Cauchy data and an improved version of the 
latter under some additional regularity hypotheses. We then use them for the 
proof of Theorem l2.4l In section 4, we derive several versions of the three spheres 
inequality (Theorems 14. 1( l4~2l and Proposition I4.5|) . which are needed to prove 
Proposition 13. II Section 5 is devoted to the proof of Proposition 13.31 which will 
use an estimate of continuation from Cauchy data, Theorem 15.21 which will be 
proven in Section 6. 
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2 The stability result. 
2.1 Notation and definitions. 

Let x e W l . We will denote by B p (x) the ball in M™ centered in x of radius 
p. We will indicate x — {x\, . . . , x n ) as x — (x',x n ) where x' — (xi . . . x n -i). 
Accordingly, B' p {x') will denote the ball of center x' and radius p in We 
will often make use of the following definition of regularity of a domain. 

Definition Let fi C 1" a bounded domain. We say r C dfl is of class C fe,Q 
with constants po, M > 0, where k is a nonnegative integer, a <G [0, 1) if, for 
any PeT there exists a rigid transformation of coordinates in which P = 
and 

nnB po (0) = {{x',x n ) G B po (0) s.t. x n > ip{x')}, 
where ip is a real valued function of class C k,a (B' po (0)) such that 

¥>(0) - 0, 
V<p(0) =0, if k > 1 
IMIc"°.«(b;, o (o)) < M po- 

When fc = 0, a = 1 we will say that T is of Lipschitz class with constants po, 
M . 

Remark We normalize all norms in such a way they are all dimensionally 
equivalent to their argument and coincide with the usual norms when p = I. 
In this setup, the norm taken in the previous definition is intended as follows: 



IMIc*.«(s; (0)) -^Poll^VIU-(s; o (o))+pS +a |^VU,s; o (o), 

1=0 

where | • | represents the a-Holdcr seminorm 

, nfc , mm \D k <p(x') - D k <p(y')\ 

\ D <P\a,B'(0) = SUp — — , 

x',w'eB io (o),xVw' W -V'\ a 
and D k <f = {D^ip}\p\ =k is the set of derivatives of order k. Similarly we set 

IMIL(O) = ^ I u 2 
Po Jn 



Po K Jn Jn ' 



Po 

The same goes for the trace norms ll u ll H ^ af2 j an d the dual norms ||'ii||H-i(n) ; 
llull _ i ,„ v and so forth. 

11 "H 7 (dQ.) 
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2.2 A priori information. 

Here we present all the a priori hypotheses we will use along the paper. 

(1) A priori information on the domain. 

We assume SI C M™ to be a bounded domain, such that 

<9f2 is connected, (2.1) 

with a sufficiently smooth boundary, i.e., 

dCi is of class C 2 ' a of constants po, M , (2.2) 

where a G (0, 1] is a real number, Mo > 0, and po > is what we shall treat as 
our dimensional parameter. In what follows v is the outer normal vector field 
to dtt. We also require that 

|0| < Mi/#, (2-3) 

where Mi > 0. 

We choose an open and connected portion T C d£l as being the accessible part 
of the boundary. We assume that there exists a point Po € T such that 

dnnB po (p ) cr. (2.4) 

(2) A priori information about the obstacles. 

We consider D C fi, which represents the obstacle we want to detect from the 
boundary measurements, on which we require that 

fi \ D is connected, (2.5) 

dD is connected. (2.6) 
We require the same regularity on D as we did for il, that is, 

dD is of class C 2 ' a with constants po, Mq. (2.7) 

In addition, we suppose that 

d(D,dQ)>p . (2.8) 

(3) A priori information about the boundary data. 

For the Dirichlet-type data g we assign on the accessible portion of the boundary 
T, we assume that 

suppg cc r. 

We prescribe the following compatibility condition (which is necessary for the 
existence of the solution, and is a consequence of the incompressibility condi- 
tion): 

gds = 0. (2.10) 



/, 

Jan 
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We shall also assume an apriori bound on the regularity of the flow, by requiring 
that for a given constant E we have 

\\g\\c^{T) < £■ (2.H) 

We also specify a bound on the oscillation of the boundary data g by requiring 
that, for a given constant F > 0, 

\\9\W{V) > F. (2.12) 

Note that (|2.11l) and (|2.12l) combined yield an a priori frequency type limitation 
of the form 

Under the above conditions on g, one can prove that there exists a constant 
c > 0, only depending on Mo, such that the following equivalence relation 
holds: 

-||5llci.= (r) < llffll H ion) - c Mlc*.«(r)- (2.13) 
2.3 The direct problem 

In this section we collect some useful results regarding the direct problem of 
finding weak solutions of (jl.ljl . We begin by an existence result (which is a 
classical result and can be found in [11], [19]) : 

Theorem 2.1. Let g <E C 1 ' Q (9fi) satisfy (gS]) and (pTTUl) , let !!cl"fl bounded 
set satisfying (|2.i p -(|2.4 |) . and let D C Q satisfy (|2.5 p -(|2.8 p . Then for every 
p, > there exists at least one solution (u,p) € H 1 (i7 \ D) x L 2 (57 \ D) of (jl.ljl . 

ll«l|Hi ( n\D) <C||5llci.-(r), (2.14) 
with C only depending on \x, Mq, M\. 

The following is a consequence of the above theorem combined with a global 
regularity result for solutions of the Navier- Stokes equations, see, for exam- 
ple, [H][VIII, Corollary 5.2]. We point out that an a priori higher degree of 
integrability of weak solutions has to be assumed. 

Theorem 2.2 (Regularity of solutions). Assume that the hypotheses of the 
previous theorem are satisfied and suppose, in addition, that (|2.1ip is satisfied. 
Let u be the weak solution to (jl.ip in f2 \ D, and suppose that u e H 1 (f2 \ D) n 
L"(fi \ D). Then for all < a < 1 we have that u € C 1 - a (fi\L») and 

IN C1 .<W) < CS, (2.15) 
where C > only depends on [i, a, Mq. 

Remark The requirement that u £ H 1 (0 \ D)riL n (Cl\D) is actually redundant 
(in the sense that it follows from u E H 1 (il \ D) ) when n < 4, due to the 
Sobolev embedding theorems. 
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The uniqueness issue for the direct problem is more subtle. Unlike the linear 
Stokes equations, here uniqueness is not guaranteed in general. In fact, several 
examples can be built even in low space dimensions (see |19j). As far as the 
inverse problem is concerned, whether or not the solution of the direct problem 
is unique is not relevant, for by formulating the inverse problem we implicitly 
select one particular solution of the direct problem to work with. To guarantee 
uniqueness, one can either a priori bound the norm of the solution (see pTj . 
Theorem VIII. 2.1), or take "not too large" boundary data g and viscosity fi, as 
stated by the following (which we will state, for simplicity, only for n < 4; see 
[T5] . Theorem 1.6, pg. 120 for a proof): 

Theorem 2.3 (Uniqueness for small data). Let Q and g satisfy the same hy- 
potheses of Theorem \2.1l and let u be a solution of (11.11) given by Theorem ] 2. 11 
Then: 

1. There exists a constant C\ = Ci(fi, 51) such that if w is another solution 
of (jl.ip and \\u\\ni (p\D) < Ci then w = u. 

2. There exists a constant C2 = Ci^i, f2, n) such that, if ||<?||c 1 '°(r) < C2, 
then u is the unique solution of 

3. There exists a constant C3 = C^{^l,n, HgHc 1 ^^)) such that if fi > C3 then 
u is the unique solution of (jl.ll) . 

2.4 The main result. 

Let n C R n , and T C dfl satisfy ([2T2> (|2l ]) . Let A C ft, for i = 1,2, satisfy 
(|2.5p - (|2.8p . and let us denote by ft; = fl \ Di. We may state the main result as 
follows. 

Theorem 2.4 (Stability). Let g £ C 1,Q (L) be the assigned boundary data, 
satisfying (j2l?l) - (j2~T2l) . Let u % G H x (fij) n L" ^) solve ([TT]) for D = A- If, 

for e > 0, we have 

p \\o-(ui,pi) ■ v - o-{u 2 ,P2) ■ ^|| H _i (r . < e, (2.16) 

then 

d n (dD 1 ,dD 2 )<p u](e), (2.17) 
where oj : (0, +00) — > R + is an increasing function satisfying, for all < t < -: 

Lj(t)<C(log\logt\)-P. (2.18) 
The constants C > and < fj < 1 only depend on fi, n, Mq, M\, £ and F . 
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3 Proof of Theorem 12.4 



In what follows, we shall repeatedly use the following notation, for various 
choices of the open set il and the parameter h > 0: 

fi/, = {x £ such that d(x, dil) > h}. 

The proof of Theorem 12 .41 relies on the following sequence of propositions. 



Proposition 3.1 (Lipschitz propagation of smallness). Let E be a bounded 
Lipschitz domain with constants Pq, Mq, satisfying (12.31) . Let u be a solution to 
the following problem: 

div a (u,p) = (it • V)u in E, 

divw =0 in E, (3.1) 
u = g on dE, 



g£C 1 > a (dE), g^Q, (3.2) 
gds = 0, (3.3) 



where g satisfies 



IdE 

IMIci.-0J3) < S, (3.4) 

\\9h*(dE) > F, (3.5) 

for given constants £ > 0, F > 0. Also suppose that there exists a point P G dE 
such that 

.9 = on dEnB Po (P). (3.6) 

Then there exists a constant s > 1, depending only on n and Mq such that, for 
every p > and for every x £ E sp , we have 

\Ju\ 2 dx >C P I \\7u\ 2 dx. (3.7) 

B p (x) JE 

Here C p > is a constant depending only on n, Mq, Mi, F, £ , po and p. The 
constant C p can be written as 



c 



p 



exp 



(3.8) 



where A, B, C > only depend on n, Mq, M%, F and £ . 



Proposition 3.2 (Lipschitz propagation of smallness up to boundary data). 
Under the hypotheses of Theorem \2.4\ for all p > 0, if x £ {Qi)( s +i)p> we have 
fori = 1,2: 

-4^ f \Vui\ 2 dx>C p \\g\\&.«<T), (3-9) 
Po Jb p (s) 

where C p is as in (13.81) (with possibly a different value of the term C), and s is 
given as in Proposition \3.1\ 
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Proposition 3.3 (Stability estimate of continuation from Cauchy data). Under 
the hypotheses of Theorem \2.J\ we have 

4=2 I IV^iI 2 < Cco(e) (3.10) 

Po JD 2 \D t 

4=2 I |Vu 2 | 2 < Cu,(e) (3.11) 
Po Jd!\d 2 

where u) is an increasing continuous function, defined on K + and satisfying 

w(t)<C(log|logt|)~ c (3.12) 

for all t < e , where C only depends on /i, n, Mq, Mi and £, and c > only 
depends on n. 

Proposition 3.4 (Improved stability estimate of continuation). Let the hy- 
potheses of Theorem \2.J\ hold. Let G be the connected component of fl± f~l &2 
containing T, and assume that dG is of Lipschitz class of constants po and Mo, 
where Mo > and < po < po ■ Then (|3.10p and (|3.11[) both hold with uj given 
by 

w(t) = C\\ogt\-i, (3.13) 

defined for t < I, where 7 > and C > only depend on fi, Mq, Mq, Mi, £ 
and^. 

Po 

Proposition 3.5. Let f2i and ^2 two bounded domains satisfying (|2.2|) . Then 
there exist two positive numbers do, po, with po < po, such that the ratios 
=j| only depend on n, Mq and a such that, if 

d u {vr u vr 2 ) < do, (3.14) 

then there exists Mq > only depending on n, Mq and a such that every 
connected component of f2i H f^2 has boundary of Lipschitz class with constants 
Po, M . 

We will give a sketch of the proofs of Propositions 13.11 and 13.21 in Section 
4, while Propositions 13.31 and 13.41 will be proven in Section 5. The proof of 
Proposition 13.51 is purely geometrical and can be found in [2] . 

Proof of Theorem \2.J\ Let us call 

d = d H {dD l ,dD 2 ). (3.15) 

Let 77 be the quantity on the right hand side of (|3.10[) and p. lip , so that 




(3.16) 
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Without loss of generality, assume that there exists a point X\ € dD\ such that 
distal, dD-i) = d. We distinguish two possible situations: 

(i) B d ( Xl ) C D 2 , 

(ii) n D 2 = 0. 

In case (i), by the regularity assumptions on dD\, we find a point x 2 G -D2 \ D± 
such that B tc i(x 2 ) C D 2 \Di, where t is small enough (for example, t = } 

suffices). Using (|3.9[) . with p = — we have 



|Vui| 2 da; > 

s p (^2) exp 



Am 



By Proposition 13.31 we have: 



Mlci.°(r)- 



(3.17) 



w(e)> 



exp 



(r) 



(3.18) 



which, once we recall (|2.12|) and solve for d, yields this estimate of log-log-log 
type stability: 

d<Cp {log[log|log e |]}~*, (3.19) 

provided e < e~ e : this is not restrictive since, for larger values of e, the thesis 
is trivial. If we call do the right hand side of (|3.19|) . we have that there exists 
eo only depending on n, Mq, Mi and F such that, if e < eo then <i < do- 
Proposition [33] then applies, so that G satisfies the hypotheses of Proposition 
13.41 This means that we may choose lu of the form (|3.13[) in (|3.18[) . obtaining 
f|3 . 10|) . Case (ii) can be treated analogously, upon substituting u\ with u 2 . □ 



4 Proof of Proposition 13.1 , 

In a recent paper by Lin, Uhlmann and Wang f |14j). the validity of the three 
spheres inequality has been extended to solutions u = {u\, . . . ,u n ) of Stokes 
systems of the form 

f Au + A(x)- Vu + B(x)u + Vp = 0, , , 

\ divw = 0, 1 } 

where A(x) is a measurable vector, B(x) is a measurable matrix, both satisfying 
appropriate bounds, and it is agreed that A ■ Vu = (A ■ Viti, . . . , A ■ V«„) (the • 
on the right hand side being the ordinary vector dot product). In what follows, 
it shall be convenient to write out the first equation in (|4.1[) component-wise 
(here and everywhere else in the paper, we use the convention of summation 
over repeated indexes): 

Am + aj-^- + hjUj + — — = 0, for i = 1, . . . , n. (4.2) 
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An analogous result when B = was obtained by Fabre and Lebeau [TU] and 
Regbaoui [18] . The result in its most general form allows A and B to have some 
type of singularity; we shall only state a simplified version: we will assume that 

\\A\\ cl , HBn) + \\B\\ c „, a{BR} <£. (4.3) 

Then we have: 

Theorem 4.1 (Three spheres inequality). Let u G ~H}{Br) be a solution to 
(|4.ip in a ball Br. Suppose that the functions A(x), B(x) are measurable and 
such that (|4.3I) holds. Then there exists a real number "d* 6 (0, e^ 1 / 2 ), depending 
only on n, such that, for all < r± < r2 < i9*f3 with r% < R we have: 

where the balls B Ti are concentric with Br and S € (0,1) and C > are con- 
stants depending only on £, n, ^ and 

We will also need to formulate a three spheres type inequality for the first 
derivatives of u: 

Corollary 4.2. Let u G H 1 (Br) be a solution to (|4.ip . and suppose that (|4.3p 
holds. Assume furthermore that B(x) = 0. XTien we ftawe that for all < n < 
r 2 < i?*r 3 w'tt r 3 < R: 




where <&* is the same as in Theorem \4-l\ 5 G (0,1) and C > are constants 
depending only on £, n, ^ and and t/ie 6a/Zs i? ri are concentric with Br. 

Remark We point out that the first equation (|1.1[) may be written in the form 
(|4.2[) by writing the nonlinear term as in (|1.2j) and calling (A(x))j = Mj and 
(B{x))ij — 0. Therefore, the three spheres inequalities (|4.4[) and ()4.5j) may be 
applied to solutions of in a domain as long as (14. 3[) holds in S with the 
aforementioned choices of A and B. 

Finally, we would like to recall the following propositions. The first is a version 
of Poincare inequality dealing with functions that vanish on an open portion 
of the boundary (see [15], or [3] for a precise evaluation of the constants in 
terms of the Poincare constant of the domain and the measure of the portion 
of the boundary of the domain where the function vanishes on) . The second is 
a Caccioppoli-type inequality for (|4.ip . which can be found in [T2"] . 

Proposition 4.3 (Poincare inequality). Let E C K." be a bounded domain with 
boundary of Lipschitz class with constants po, Mq and satisfying (|2.3p . Then 
for every u G H 1 (E) such that 

u = on dEr\B po (P), 
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where P is some point in dE, we have 

< C/Oo||Vu|| L 2 (E) , (4.6) 
where C is a positive constant only depending on Mq and M\. 

Proposition 4.4 (Caccioppoli inequality). Let u E H 1 (Bn) be a solution of 
(|4.ip in Bji, and suppose that (|4.3|) holds. Then there exists C > depending 
only on n,/j, and £\ such that, for every r with < r < R we have 

|Vu| 2 < — / | U | 2 . (4.7) 



These two facts can be used to prove Corollary 14.2 



Proof of Corollary \4-2\ We notice that if u is a solution to (|4.1|) and B(x) = , 
then u — ue is a solution as well, where ue denotes the average of u in E. Since 
the classical Poincare- Wirtinger inequality (see [TS] ) holds for u — ue, we apply 
it together with the Caccioppoli inequality and the three spheres inequalities 
(|4~4l) . (|43]l to obtain: 



( r3 _ r2 )2 



B,, 



and Ci, C2, C3 only depend on £, n, ^ and £|. □ 

We can now show that Proposition 13.21 follows from Proposition 13. II 
Proof of Provosition HOI From Proposition 13.11 we know that 



V Ml | 2 > C p / _|Vu,| 2 , 

B p (x) Jn\Di 



where C p is given in (|3.8I) . We have, using Poincare inequality (I4.6[) and the 
trace theorem, 

/ a _|v^| 2 > CtfT'WM^am > WMUw (4 ' 8) 

Applying the above estimate to (|3.7|) and using (|2 . 1 3|) will prove our statement. 

□ 

Finally, we will need a three spheres inequality for functions that can be 
written as differences of solutions of the Navier-Stokes equations (due to non- 
linearity, this does not follow the previous remarks). 
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Proposition 4.5. Let u\ and u 2 be solutions of 



div a{u t ,pi) ={ Ui -V)ui in E, . . 

div Ui — in E, 

fori = 1,2. Suppose that ||mi||c 1 - q (£;) + II ^2 1| c°>" ( _E) < £• LetBn(x)cE. Then 
there exists a real number € (0,e -1 / 2 ), depending only on n, such that, for 
all < n < ^2 < i9*r3 wif/i r% < R we have, calling w = U\ — u 2 : 

8 / r \ 1-5 



\w\ 2 <c( \w\ 2 ) ( / | W | 2 ) " (4.10) 



where the balls B ri are centered in x, and S £ (0, 1) and C > are constants 
depending only on £ , n, ^ and 

Proof. In view of Theorem (|4.ip (and the subsequent remarks), it is enough to 
show that w can be written as a solution to a system of the form (|4.1[) . This is 
readily done, by subtracting from each other (|4.9[) for i — 1,2. We may write, 
in components for j = 1, . . . , n: 

, , duij d(ui)j d(pi — P2) 

Calling (A(x))i = {-u 2 )i and (S(a;))y = we have that holds in E, 

so the hypotheses of Theorem 14 . 1 1 hold for w. □ 

Remark We observe that, since the identically zero function solves (|4.9[) in 
E, we can also apply the three spheres inequality to each Ui separately (as we 
already pointed out in the previous remark). 

The proof of Proposition 13. II is now a consequence of the work done so far. 



Proof of Proposition \3.1\ The proof is based upon the validity of the three 
spheres inequality for 11.11 Since we have just proved Proposition I4.5( we re- 
fer to [7], [17] and [16] for a detailed proof, which applies here with only slight 
modifications, as it only requires (|4.4[) and (|4.5|) and some geometric construc- 
tions which only exploit the regularity of dE and thus apply unchanged to this 
situation. □ 



5 Stability of continuation from Cauchy data. 

We will need the following lemma, proved in [2]: 

Lemma 5.1 (Regularized domains). Let be a domain satisfying (|2.2p and 

(|2.3I) . and let Di, for i = 1,2, be two connected open subsets of fl satisfying 

(|2.7p . (|2.8[) . Then there exist a family of regularized domains C tt, for 
< h < apo, with C 1 boundary of constants po, Mq and such that 

Di C C D'l 2 ifQ<h x < h 2 ; (5.1) 
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j h < dist(x, dDi) < j x h for all x £ dD%; (5.2) 

measp? \ A) < l2M lP 2 h; (5.3) 

meaSn-x^Df) < j 3 M lP 2 ; (5.4) 
and /or every x £ dD\ there exists y £ dDi such that 

h a 

\y-x\ = dist(x,9A)j K 33 ) - Ky)l < 74— ; (5.5) 

Po 

where by v(x) we mean the outer unit normal to dD^, v{y) is the outer unit 
normal to Di, and the constants a, ^j, j — . . .4 and the ratios j^, ^ only 
depend on Mq and a. 

The proof the stability estimate of continuation from the Cauchy data heav- 
ily relies upon the upcoming result, which deals with the estimation of the 
stability of the stationary Navier-Stokes equations with homogeneous Cauchy 
data, the proof of which, in turn, is based upon an extension argument. We 
pospone the proof to the next section. Here we will set up the notations and 
state the theorem. Let us consider a bounded domain E C M" satisfying hy- 
potheses (|2.2p and (|2.3I) . and take T C dE a connected open portion of the 
boundary of class C 2 ' Q with constants po, Mq. Let Po G T such that (|2.4I) holds. 
By definition, after a suitable change of coordinates we have that Po = and 

EnB Po (0) = {(x',x n ) £ E s.t.x n > ip(x')} c E, (5.6) 

where ip is a C 2,a (B' (0)) function satisfying 



*»(0) = 0, 
|Vv>(0)| = 0, 
Mlc 2 .°(B' (0)) < M po- 



Define 



Poo 



V 1 + M l (5.7) 

r = {(x',x n ) £ r s.t. \x'\ < poo, x n = <p(x')}. 

In what follows we shall consider (ttj,Pi), for i = 1,2, which are solutions to the 
following problems: 

div a(ui,pi) =(ui-V)ui in E, 
divu t =0 in E, 
Ui = g on 1 , 
o-(ui,Pi)-v =ipi on T, 

where g satisfies (|3.2p - (|3.5l) . ipi £ C 1,a (r) and we use the same notations as 
we did in Section 4 (which are also to be understood in what follows). Define 
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w = u\ — u 2 , q = pi — p 2l these will solve a system of the following form: 

div cr(w, q) —A-Vw + Bw in E, 

div w = in E, . 

w = on r, ^ > 

<j(w, q) ■ v — ipo on r, 



where A and B were explicitated in (|4.2I) and -0o = "01 — ^2- We have the 
following estimate for a solution of systems of the form (|5.9p . the proof of which 
is delayed to the next section: 

Theorem 5.2. Let Abe a vector ofC 0,a (E) of constants po, Mo, let B a matrix 
of class C 0,a (E) of constants po, Mo, satisfying 

\\A\\co.- + \\B\\co,«=e, (5.10) 

and let (w,q) be a solution of class C l a (E) x C°(E) to the problem: 

div a(w,q) =A-Vw + Bw in E, 

div w = in E, , . 

U) =0 on r, l& j 

cr(w, g) • ^ = "0 on r, 

where tp £ H _ 2(L) and T C BE is of class C 1,a . Then there exists p, only 
depending on Mo, a and £, such that, letting P* — Pq + jv where v is the outer 
normal field to dE, we have: 

IHIl^hb^ (p-)) < 4lhllL2 ( T B) (po||^ll H -^ r J r , (5.12) 

where r only depends on a and Mo and C > also depends on £ . 

Proof of Provosition VS . 3[ Let 9 = min{a, ^ 2 ~ I0 (i+m' 1 ) } wnere °j 70j 7i are 
the constants depending only on Mo and a introduced in Lemma 15.11 then 
let p = 9po and fix p < ~p. Let D^, D 2 be the regularized domain associated 
with Di, D 2 respectively, built according to Lemma IBTTl Let G be the connected 
component of 0\ (D\ U D 2 ) containing 9fi, and G p be the connected component 
of O \ (D^ U D 2 ) which contains <9f2. We have that 

D 2 \ThcQi\Gc ((£>? \5I)\G)u((n\G")\I>f) 

and 

fl((n\G")\i^) = rfur§, (5.13) 

where L£ = <9L>£ n dG p and rf C dDf . Then 



Vwi| 2 < / |Vwi| 2 < / |Vui| 2 + / |Vmi |2 



_ _ 1 

D 2 \D 1 Jn x \G J(DP\Di)\G J(Q\Gp)\D? 



(5.14) 
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The first term can be estimated directly, using (I2.15[) and (|5.3[) we have 

/ _ _|v Ml | 2 <c Po "- 2 || 5 ||^, Q(r) ^- (.vi.r, 

J(D?\D^\G PO 



where C only depends on the M , M X) a and p. We call Cl(p) = (fi \ G p ) \ Z?^. 
We use the first equation in multiply it by u\ and integrate over VL(p) to 

derive the following identity: 







/ [iui ■ Aui — ui(ui ■ V)ui — u\ ■ Vpi = 
Jn{ P ) 

\iu\ ■ Aui — -div (ui|ui| 2 ) — u\ ■ Vpi — 
n( P ) 2 

= / p(V«i • v)u\ - -Mi • z/|wi| 2 - (ui ■ v)p - I |Vui| 2 , 
JdQ(p) 1 Jn{p) 

which yields, recalling (|5.13|) . 

/i / |Vui| 2 = / ^(V-ui • v)m - ^ui • v\m\ 2 - (ui • v)p + 

Jn[p) y \ (5.i6) 

We start by estimating the first integral on the right hand side of (|5 . 16[) . If 
x € r p , by Theorem 15. 1[ we find y G dD\ such that \y — x\ = d(x, dDi) < jip; 
since u\(y) — 0, by Lemma 12.21 we have 

\ Ul (x)\ = \ Ul (x) - Ul (y)\ < C^\\g\\ c i, a(r) . (5.17) 

Po 

On the other hand, if a; £ T^, there exists y S D 2 such that \y — x\ — d(x, dD 2 ) < 
7ip. Again, since u 2 (y) = 0, we have 

\ui(x)\ < \u x {x) - ui{y)\ + \ui(y) - u 2 (y)\ 
< C (A Wol .. (r)+agfKiH) , (5-18) 

where w — ui — u 2 . Combining (|5.17[) . (15.181) and (|5.16[) and recalling (|2.15p 
and (15.41) we have: 

/ |Vui| 2 < — Po~ 2 (ll5llci.°(r)+l+M)(lMlci>=(r)— +Nld.<*(r) max Ju; 
Jd 2 \d 1 P k Po dcp\an 

(5.19) 

We now need to estimate maKg GP \g n \ w\. We will do so by means of Proposition 
I4~5l Take x e 9ff \ <9f> and define 

o* = min / . Po — ^-^l, (5.20) 
P \ 16(1 + Af 2 ) ' 16 J ' v ; 



16 



min{p 00 ,p} 

x Q = P Q — v, (5.21) 



where v is the outer normal to dfl at the point Pq. By construction io € O f 



There exists an arc 7 : [0, 1] H> G p H flp^ such that 7(0) = xq, 7(1) = a; and 



7([0, 1]) C G p n . Let us define 

2 

3 1 

p 3 = min{7optf*,p}, p 2 = ^3, Pi = -^p3, (5.22) 

where I?* is the constant given in Theorem 14.11 We pick a sequence of S + 1 
times ti and points Xi = j(U), i = . . . S, defined by the following construction. 
Call t = 0, then: 

U = max{i e (0, 1] s.t. | 7 (f) - ^ I = y }, if |ar< - z| > y , 
otherwise, i = S, 
and stop the process. The number of spheres is bounded by 

where C only depends on a, Mq, Mi and E. The balls Bp3_(xi) are pairwise 
disjoint, the distance between two consecutive centers is given by 

|x i+ i - Xi\ = y , i = 0...S-l, |xs-a;|<y. 

We iterate the three spheres inequality (|4.10[) on a chain of spheres with radii 
pi, p 2 and p 3 , this leads us to 

\w\ 2 dx <C( \w\ 2 dx) ( \w\ 2 dx) (5.23) 

where < S < 1 and C > only depend on Mq, a and £. From our choice of 

p and 1?*, it follows that B pi (x ) C B p *(x ) cGfl P>3pi (P*), where we follow 

4 

the notations from Theorem 15.21 We may apply Theorem 15.21 to w; thus, using 
(|5T2|) . dHHJ, ([2~T6]) and ([2~TT|) on ([5^3)) we have: 

\w\ 2 dx < Cp™- 2 e 2TsS . (5.24) 

B P2 (x) 

The following interpolation inequality holds for all functions v defined on the 
ball B t (x) C W l : 



(5.25) 
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We apply it to w in B P2 (x), using (|5.24|) and (|2.15p we obtain 

IHIl~ (Sp2 (,))<c(^)V S , (5.26) 

where 7 = ^q^. Finally, from ()5.26|) and (|5.19|) . and recalling (|2.11j) we get: 

/ |V Ul | 2 <Cpr 2 (^+(^) 5 e^ S ), (5.27) 

with C only depending on Mo, a, p and £. Let us now choose p depending 
upon e, of the form 



2Slog|S| 
log I log eT| 



We have that p is defined and increasing in the interval (0, e 1 ). Call ( the 
number such that p(() = min{p, />}, and let C = m in{C? exp(— 7 2 )}. Since the 
thesis is trivial for larger values of e, it is not restrictive to prove it only in the 
smaller interval (0,£). We are able to apply (|5.27p to (I5.14|) with p = p(e) for 
e G (0, £) to obtain 



/ |V Ul | 2 < Cp^log I log £P, 

JD 2 \D! 

and since e < exp(— r y 2 ) it is elementary to prove that 



(5.28) 



log I log e 0, 1 > -log|loge|, 



so that (|5.28|) finally reads 

/ |V«i| a < CprM*), 
Jd 2 \d 1 

with u{t) = log I logtj n defined for all < t < e -1 , and C depends on M , Mi, 
a, p and £ . □ 



Proof of Proposition \3.4] We will prove the thesis for u\, the case 112 being 
completely analogous. First of all, we observe that 

/ p|Vui| 2 < / p|Vui| 2 = / ^(Vu r ^)ui-pi(u r ^)-iu r ^|ui| 2 , 

JD 2 \D 1 Jnx\G Jd(Ui\G) 1 

(5.29) 

and that 

«9(fii \ G) C 9-Di U (<9L> 2 n 8G). 

If we recall the no-slip condition, apply to (j5 -29[) computations similar to those 
in (JEHJ), (jS35]) . we get 



/ |Vui| 2 < -p£ 2 (||5llci.= ( r) + 1 + M) max |w| 
Jd 2 \d 1 M dD 2 ndG 



(5.30) 
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where again w = u\ — U2, and C only depends on a, Mq and M\. Take a point 
z € 5G. To evaluate max \w\ on <9D 2 H <9G, we start by choosing a point z 6 9G 
and estimating ||Vu|| on a ball centered in z, in terms of ||Vtt|| evaluated on 
a ball centered in xq- We will do so by applying iteratevely the three spheres 
inequality, twice. By exploiting the regularity assumptions on dG, we find a cone 
centered in z, which we denote by C(z, £, $0) (where d = arctan is half the 
aperture of the cone and £ € K n is a unit vector representing the direction of 
the cone), such that C(z,^,'3o) H B Po (z) C G. It can be shown (0 Proposition 
5.5]) that G p is connected for p < ^y 2 - with /i only depending on M . We now 
claim (without proof, see [7] and [TB] for the detailed constructions in the same 
context) that we may build Ai > and Q\ > 0, such that, if we define 

wi = z + Ai£, 

pi = -d*h \i sintfi. 

where < d* < 1 was introduced in Theorem 14.11 then the following claims 

hold: B P1 (wi) C C(z, £, J?i) n (z) and B Vl (iu x ) C G(z, £, i? ) n 5^ (2) C G. 

$ * 

Furthermore ^£f- < P*, hence J54 P1 (x ) C G, where p* and xo were defined by 
v w* 

(|5.20[) and (|5.2ip respectively. It follows that w\, xq € Gi? x , which is connected 
by construction. Iterating the three spheres inequality (|4.10l) (mimicking the 
construction made in the previous proof) 

/ \w\ 2 dx <C( \w\ 2 dx) ( / \w\ 2 dx) (5.31) 

JB n ( Wl ) Wg ' W.B P1 (xo) ' 

where < S < 1 and G > 1 depend only on n and £, and S < We apply 

Theorem 15.21 in the same fashion as the previous proof, which leads to 

/ h 2 < c P y?, 

Jb p1 ( Wi ) 

where < (3 < 1 and C > only depend on a, Mo, £ and So far the 
estimate we have is only on a ball centered in wi , we need to approach z € 9G 
using a sequence of balls, all contained in C(z, £, #1), by suitably shrinking their 
radii. Take 

1 — smi?i 
X ~ 1 + sin^i 

and define, for k > 2, 

Afc = xA/t-i, 
Pk = XPk-l, 

w k = z + A fc £. 

With these choices, A fe = Ax fe_1 Ai, p k = X* -1 pi and B Pk+1 (w k +i) C B 3pk (w k ), 
B^ pk (w k ) C C(z, £, i? ) n B po (z) C G. Denote by 

d(k) = \w k - z\ - p k , 
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we also have 

d(k) = x*- x d(l), 

with 

d(l) = Ai(l-?9* sini?i). 

Now take any p < and let fc = the smallest integer such that d(k) < p, 
explicitly 

l^< fc(p) _ 1 <^i*l 

logx logx 
We iterate the three spheres inequality (14.5[) over the chain of balls centered in 



< Hp) - 1 < , — + 1. (5.32) 



Wj and radii pj, 3pj, 4pj, for j = 1, . . . , fc(p) — 1, which yields 

/ M 2 < C P ™e 2 ^ (P, ~\ (5.33) 

with C only depending on a, Mo, £ and Using the interpolation inequality 
(pT2"5j) and ([2TT5]) we obtain 

e /3i<5 fc(p >- 1 

lkll^(B Pfc(p) (^ (p) )) < P y!(tM _ ir (5-34) 

where A = ^ depends only on a, M , Mi, £ and From ([5~3"Ij) ) and (|2~T5|) 
we obtain 

Finally, call 

p(6)=rf(l)|l0g6^r S , 

with 

I log X I 



B 



2 log |£| 



and let £ = exp(— f3 1 1 ). We have that p(e) is monotone increasing in the interval 
< e < (, and p(C) = d(l), so p(e) < c?(l) there. By choosing p — p(e) from 
(jQ5j) and (jOO)) we obtain 



/ |V Ul | 2 <C^ 2 |loger s , 



(5.36) 



where C only depends on fi, a, Mo, £ and □ 
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6 Proof of Theorem 15.2 



This section is entirely devoted to the proof of Theorem 15.21 Let us start by 
introducing some notations. We define 



Q(Po) = B' poo (0) x 

and 



m qP 1 MqpI 



T Q = dEnQ(P ). (6.1) 

Finally, let us call E~ = Q(P ) \ E and E = E U E~ U T . We start the proof 
by choosing a vector A and a matrix B such that A — A, B — B in E, and 

WMc^(E) + l|5|U«(£) ^ c i £ > ( 6 - 2 ) 

where Ci > only depends on a and Afo- Our aim is to build an extension w 
of w such that it satisfies the extended problem 



div a(w,q) — A ■ Vw + Bw + $ in i?, 
div w = in E, 



(6.3) 



where $ <E H 1 (E) is such that 

||$|| H _ I(ii) <qMI H _ i(r) . (6.4) 

Define u>, q as follows: 

u> in £7, 
in E~. 

q in E, 
q = \ in 

We have that w € H x (i?) and that div w = 0, in the weak sense in E. 

In order to write a system (|6.3p for we take any i> G Hj(i?) and 

consider 



er(ui, (?) ■ Vu = / cr(w, q) ■ Vv + / a(w, q) ■ Vu. (6-5) 
By the divergence theorem on the first term we obtain 

/ a(w, q) ■ Vv = - / (A - Vic) • Vv - / Bw ■ Vv + / ip -v. (6.6) 
J E J E J e Jr 

Define $(w) = J T ip ■ v for all v £ Hq(E). Using (|6.6p and the trace theorem: 

|*(»)| < U\\n- kr) \\v\\ Hkr) < CA,|MI H - i(r) IMlH W , (6.7) 
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therefore w satisfies (|6.3[) in the weak sense, with $ such that (|6.4[) holds. We 
now want to apply the three spheres inequality to the inhomogeneous system 
(|6.3I) . In order to do so, we need to establish local well posedness for the Cauchy 
problem for the linearized Navier Stokes equations. We claim the following: 
there exists p such that the problem 

div a(w*,q*) = A ■ Vw* + Bw* + $ in Bp, 

divw* =0 in Bp, (6.8) 
w* = on dBp. 

admits a weak solution w* in the ball B-r~ such that 

rno 

\\™*\\Kl(B ? )<Cpo\m\H-HE-), (6.9) 

This can be shown by projecting (|6.8I) on the space of divergence free function, 
so that it becomes a pressure free second order partial differential equation in 
w* only, with the principal part being the laplacian operator, and the lower 
order terms are continuous and bounded by a constant depending only on a, 
M and 8. Therefore the existence of a solution and the well posedness of (|6.8p 
is shown by a standard coercivity argument for sufficiently small radii. Using 
linearity we write w = wo + w* , with w* € Hq(Bp) such that (w*,q*) solves 
(|6.8p and wq solves 

div a(w Q , q ) + A ■ Vw + Bw = in B pi , , g ^ 

div wo =0 on B P1 . 

Using interior regularity of solutions for elliptic systems we get 

lko|k°°(Bt(aO) <t~*\\wo\\L?(B t (x))- (6-11) 

We will thus need to estimate || w ||l 2 (b(x)) on a ball near the boundary. We 
may apply Theorem 14. II to Wq, thus, calling and r 3 = p, r\ — r 2 = we 
have (understanding that all balls are centered in Pq) 

\\wo\\v(Br 3 ) < C\\wo\\l2 {Bri) \\wo\\l^ Br3) , (6.12) 

with C > only depending on n and £. Let us call r\ = po\\ip\\ _i ■ By the 

h 2(r) 

triangle inequality we have that 

||^o||L2(B r ) < ||w||L2(S r ) + ||«>*||L3(B r ) < ||w|| L 2 (Bi . ) + Cf], (6.13) 

for r = r\ , r3 ; furthermore, 

IN| L 2(B r2 ) < |ko||L 2 (s r2 ) + l|w*|| L 2(B r2 ) < l|wo||L2(s r2 ) + C77, (6.14) 

where (in both cases) C > only depends on n and £. Putting together (|6.12p . 
([633)) . l|6"Tljl we get 

IH| L 2(B r2 ) ^ IHIi^B^nB) < 
<Ct7 + C(|H| L2(Sn) + C77) r (||«;|| Ija(Bt3nB) + < (6.15) 

<C(ry + rTfr + WwU^e)) 1 --) < C^\\w\\^ Ey 
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